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CrN is observed to undergo a paramagnetic to antiferro- 
magnetic transition with Neel temperature TV ~ 280 K, 
accompanied by a shear distortion from cubic NaCl-type 
to orthorhombic Pnma structure. Our first-principles plane 
wave calculations, based on ultrasoft pseudopotentials, con- 
firm that the distorted antiferromagnetic phase with spin con- 
figuration arranged in double ferromagnetic sheets along the 
[110] direction is energetically much more stable than the 
paramagnetic phase, and also slightly favored over the ferro- 
magnetic and other examined antiferromagnetic phases. The 
energy gain from polarization is much larger than that from 
distortion; nevertheless, the distortion is decisive in resolv- 
ing the competition among the antiferromagnetic phases: the 
anisotropy in the (100) plane arising from the ferromagnetic 
double-sheet magnetic order allows a small but important 
energy gain upon the cubic-to-orthorhombic transition. Al- 
though antiferromagnetic order leads to a large depletion of 
states around Fermi level, it does not open a gap. The system 
is metallic with occupied hole and electron pockets of Fermi 
surface containing ~ 0.025 carriers of each sign per formula 
unit. The simultaneous occurrence of structural distortion 
and antiferromagnetic order, as well as the competition be- 
tween antiferromagnetism and ferromagnetism, is analyzed. 
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I. INTRODUCTION 

3<i transition metal mono-nitrides (TMN) are inter- 
esting technologically due to a rich variety of proper- 
ties, ranging from high-T c superconductivity (observed 
in TiN and VN), to, very hard surface coatings, to het- 
erogeneous catalysisa. Although some characteristics are 
common to the whole series of compounds such as hav- 
ing high melting points and chemical hardnesses, this 
class also exhibits interesting sanation in some prop- 
erties. ScN is semiconductingcl, as CrN is sometimes 
reportedB to be, while TiN and VN are superconductors!! 
Regarding magnetic properties, ScN, TiN and VN are 
paramagnetsB, CoN has been only recently synthesized^ 
and is also found paramagnetic. For CrN, strong ex- 
perimental evidence supports antifexromagnetism below 
a Neel temperature T~273-286 KHijIj, and repent mea- 
surements suggest for FeN antiferromagnetisn£3 as well. 

Even for the structural properties there is no unifor- 
mity. Most of compounds have rocksalt structure, but 
there are exceptions: CoN has ZnS-type structures, and 
FeN has been observed in both ZnS- and NaCl-type 
phases (depending on the growth conditions). CrN seems 



to be very peculiar in that it is the only one of the series 
that does not have a ground state with cubic symmetry. 
Its transition from paramagnetism to antiferromagnetism 
at Tn is accompanied by a structural distortion from cu- 
bic to orthorhombic (Pnma) symmetry^!. 

The body of information on 3d TMN properties is 
rather fragmentary (and sometime contradictory), and 
much vital data are unavailable. Most of these com- 
pounds are difficult to synthesize, and the strong ten- 
dency of nitrogen to form dimers. leads to difficulty in 
growing defect-free TMN crystals^ (a problem also com- 
mon to other nitride compounds such as III-V semicon- 
ductors). Thus far, nothing is known about MnN, NiN, 
CuN and ZnN, and. only recently it has become possible 
to synthesize FeNtL3 and CoNR 

Due to the difficulty in preparing stoichiometric CrN, 
and to a strong sensitivity of its properties on the con- 
centration of defects, reports on the properties of CrN 
inJjhe literature are not entirely consistent. Browne et 
al.a have studied the effects of variation of x in CrN^ 
and produced samples with x up to 0.997, taking care 
to keep them oxygen-free. Only at this high N compo- 
sition is the antiferromagnetic transition sharp, and it 
occurs at Tat=286 K. The transformation is first order, 
with hysterisis of 2-3 K as measured by linear expansiv- 
ity, susceptibility, and resistivity. At lower values of x 
(down to 0.980) the transition temperature decreases to 
273 K and the transition width increases. 

At x=0.997, the properties reported by Browne et 
al. are as follows. The resistivity slope dp/dT is posi- 
tive (metallic-like) on both sides of the transition, and 
p decreases by 30-35% across the transition, in opposi- 
tion to the semiconducting behavior reported by Herle 
et al.a Due to the polycrystalline nature, and specific 
density only 70% of the bulk, the value of the resis- 
tivity ~ 2 — imTtcm is only an upper bound. The 
susceptibility slope changes from a small negative value 
dlogx/dT ~ —4 x 10~ 4 KV 1 above T^ to a positive value 
of similar magnitude below Tjy. The decrease in x is ~ 
35%. These data are most consistent with a metal-to- 
metal transition at Tjy in which a substantial fraction 
of the Fermi surface is gapped by structural distortion 
and/or the onset of magnetic order. Tsuchiya et al.tl 
have produced thin films of CrN^ with 1.0 < x < 1.2. 
These films, like the samples with x < 1 of Browne et 
al., have Tat < 280 K. None of the transitions are as 
sharp as for the x = 0.997 samples. It is clear that prop- 
erties are sensitive to N overstoichiometry as well as to 
understoichiometry. 

Not surprisingly, the theoretical literature is also not 
very extensive. Some work on the 3d TMNs in the cu- 
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bic paramagnetic phase has been reportedt2rt-3, but very 
little is reported for compounds that undergo reconstruc- 
tions accompanied by, and perhaps.-GLriven by, .magnetic 
ordering, as may be the case in CrNlij and FeNliZI. Inves- 
tigation of compounds that are so close in chemical com- 
position but with so different physical properties are im- 
portant in providing insight into basic mechanisms such 
as magnetic ordering, formation of local moments, and 
relations of this phenomena with basic features such as 
structure and density of states. Ab initio calculations, 
which can separate the effects of magnetism and struc- 
tural distortions, are very effective in such cases. In this 
paper we focus on CrN and investigate the occurrence 
of magnetic ordering for the cubic structure and for the 
(observed) orthorhombic structure. 

A basic characteristic of CrN in the cubic paramag- 
netic phase, which makes it easy to understand the in- 
stability towards ferromagnetism, is the very high den- 
sity of statesEJ at the Fermi energy (E^). Less easy to 
understand is the competition between AFM and FM or- 
dering, and the relation between spin polarization, mag- 
netic order, and structural distortion. Other cases are 
well known where high DOS drives systems to structural 
(e.g. VO2, from tetragonal to monoclinic) or spin polar- 
ized (Cr02, from PM to half-metallic FM) transitions. In 
fact, the isovalent compound MoN shows an interesting 
parallel. Due to the peak at E^ in the density of states in 
the cubic phase, it was predicted to be a high tempera- 
ture superconductorlia (compared to conventional supex=. 
conductors). It was subsequently shown by calculations 
that it is structurally unstable in the cubic structure, 
which explains the fact that experimentalists had been 
unable to grow stoichiometric samples. In the case of 
CrN the cubic structure is stable at high temperature, 
but at Tpf there is the simultaneous occurrence of both 
structural transformation and magnetic order. Giving 
insight into possible relations between these transitions, 
and identifying the driving force for the transition, is the 
main motivation of the present work. 

The paper is organized as follows: in Section |l| a re- 
view of the CrN crystal structure is given, together with 
a brief description of our calculational procedure. In Sec- 
tion III we show structural and energetic results for CrN 



in different magnetic phases. Finally, Section IV is dedi- 



cated to an analysis of the electronic properties. Results 
for density of states, band energies, and some Fermi sur- 
face properties are presented and discussed. 



II. CRYSTAL STRUCTURE AND METHOD OF 
CALCULATION 

In Figure [l] a tridimensional perspective (top panel) 
and a top view down onto the (001) plane (bottom panel) 
of the CrN structure is shown. The x, y and z axes 
of the unit cell for the AFM phase are oriented along 
the [1T0], [110] and [001] directions of the cubic NaCl 



structure, respectively. The shear distortion consists of 
a contraction (not visible in the Figure) of the angle a 
from the cubic (a— 90°) to a ~2° smaller observed value. 
The square symmetry in the (001) plane is broken so that 
the coordination of the atoms reduces to 2 in the (001) 
plane. Orthorhombic and cubic structures are connected 
by the relations a = 2aosin(a/2), b = aocos(a/2) and 
c = a . 

From experiments, two possible antiferromagnetic spin 
alignments have been proposedotl. We consider here prL. 
marily the arrangement determined by Corliss et al.u 
and depicted by arrows on Cr sites in Figure [j]: dou- 
ble ferromagnetic layers parallel to the [110] direction 
alternate in such a way to form a four-layer antiferro- 
magnetic orthorhombic unit cell (4 formula units and 8 
atoms per cell), with non vanishing magnetization field 
planarly averaged along [110]. As a result of this order- 
ing, six of the 12 (cubic) Cr neighbors have spin parallel 
and six spin antiparallel. Of the six (cubic) next near- 
est neighbors, four have antiparallel spin. Ref. ^ reports 
for the magnetic moments m = 2.36 hb per Cr atom. 
Another experiment!] gives a much larger m = 3.17 hb 
that seems to be much less reasonable. Indeed, the first 
value is consistent with a Cr to N charge transfer of ~ 2 
electrons, while the second would imply an almost com- 
pletely ionized configuration Cr + ( 3-(5 )N~( 3 ~' 5 ) that is not 
to be expected (charge transfer would be larger than in 
more ionic compounds such as transition- metal oxides). 

Another motivation for this work is technical and 
centers on our method of calculation. These local- 
spin-density calculations are performed within a plane 
wave and ultrasoft pseudopotentiacj (USPP) framework. 
Here we want to furnish a sound assessment of accuracy 
and efficiency of the use of pseudopotentials for the study 
of materials historically considered impossible, or at best 
difficult, to be treated within a plane wave pseudopo- 
tential method. CrN is a difficult test, being the com- 
bination of a 3d transition metal and a first-row atom, 
with high density of states at and a reconstruction 
accompanied by antiferromagnetism. 

To our knowledge this is the first plane-wave study 
of AFM pCrN (only one augmented-spherical-wave 
calculationllj is present in literature). The use of pseu- 
dopotentials ta_study magnetic compounds itself is only 
two years QldE3 Use of USPP for magnetic materials 
is very newESta and presents the essential feature: con- 
structing smooth (but well transferable) USPP allows us 
to obtain converged results (within 1 mRy for the total 
energy) at energy cut-off of 29 Ry (^2600 plane waves per 
CrN unit), whereas ordinary norm-conserving pseudopo- 
tentials would require at least the twice as large a cutoff 
value, hence four times as many plane waves. The chosen 
core radii are 2.20 a.u., 2.45 a.u. and 2.45 a.u. for s, p 
and d projectors of Cr pseudopotential, respectively; 1.30 
a.u. for both s and p components of N pseudopotential. 
Two projectors for each angular channel (corresponding 
to energies ^3 eV apart) were used. The (highest in en- 
ergy) p component was employed as local part. 



2 



The non-linear core correction (NLCC)c3 was used to 
ensure the best degree of transferability. As pointed out 
in Ref. [25l the NLCC is particularly important for spin- 
polarized situations. Indeed, pseudopotentials are gen- 
erated in a non-spin-polarized reference configuration. 
Such a procedure is justified only to the extent that the 
pseudopotential is perfectly descreened by the valence 
dependence. In case the overlap between core and va- 
lence charge is large, as it happens for transition metals, 
the non-linearity of the exchange-correlation potential in- 
duces incomplete descreening, and thereby a pseudopo- 
tential dependence on the valence configuration. Also, 
the LSD exchange-correlation potential is even more non- 
linear than the LDA onec3, and further affects the pseu- 
dopotential accuracy. The NLCC set the pseudopotential 
free from its exchange-correlation dependence, thus re- 
covering its transferability. Finally, no use of relativistic 
correction was made. 

The performance of Cr and N pseudopotentials has 
been individually tested. In particular, with the former 
we obtain an accurate description of bec bulk Cr (e.g. a 
calculated lattice constant within 0.5 % the experimental 
measure), whereas the latter has been estensively em- 
ployed in a long investigation of semiconducting III-V 
nitrides compounds (see ref. ^6| and references therein). 

The exchange-correlation potential was that of Cepeft- 
ley and Alder, as parametrized by Perdew and ZungerEZI. 
We sampled the irreducible zone (IBZ) with special k- 
points (~ 30 for self-consistent calculations, from 150 up 
to 530 k points to evaluate the densities of states). 

III. STRUCTURAL PROPERTIES 

In Table || our results for energies and magnetic mo- 
ments of PM, FM and AFM phases of CrN are shown. 
For a fair comparison between theory and experiment we 
investigated, together with the AFM double-sheet order- 
ing found by Corliss et al (here indicated as AFM? 110 j), 
other two kinds of AFM orderings: the single-sheet 
AFM^ 110 j phase, similar to the AFMj? 110 j but with single 
ferromagnetic sheets having alternated spin along [110] 
(a tetragonal structure with 4 atoms per cell), and the 
AFMnii] phase, consisting on single ferromagnetic sheets 
parallel to the (111) plane which alternate the sign of the 
spin along [111] (i.e. a trigonal structure common in 3d 
transition metal monoxides with 6 Cr and 6 N layers in 
the cell). For each phase, we present results obtained in 
both cubic and distorted structures. For the cubic struc- 
tures we consider both experimental and theoretical lat- 
tice constants (i.e. the ao obtained by energy minimiza- 
tion of the respective phases). To add the further effect 
of the planar distortion, we extract the value of a from 
the experimental structure given in Ref. ^ (a=88.23°). 

A number of results for the P-MpCubic (NaCl) phase 
of CrN are present in litetatureE^ESO. In particular, 
Shimizu, Shirai, and SuzukitH performed lattice constant 



calculations for the whole series of 3d TMN's. For CrN 
they report a,Q M —7 .56 a.u. which underestimates the ex- 
perimental value aQ Xp —7.81 a.u. for the AFM^ 110 j phase 

(a 10% smaller volume). Our calculation gives M =7.73 
a.u., a value 1% smaller than experiment. The discrep- 
ance between our a^ M value and that of Ref. |l7] can be 
ascribed to the different form of the exchange-correlation 
used and, moreover to the fact that, differently from Ref. 
[L~7| , we perform a non-relativistic calculation. 

For a true comparison with the experiments, however, 
it is necessary to calculate the lattice constant for the 
cubic AFMj\ 10] structure. We find afi FM =7.876 a.u., 
less then 1% above experiment. Thus, according to our 
calculations, polarization (from PM to AFM^ 110 j order) 
produces an expansion of 2% of the lattice constant. 

Looking at Table |, we see that the AFMj? 110 ] distorted 
phase is favoured overall. The theoretical determina- 
tion of ao produces, with respect to the experimental 
structure, an energy gain (E^ x -Ej") of 0.02 eV/formula 
unit; anyway, in general, results obtained for theoretical 
and experimental lattice constants are close each other to 
within 0.01 eV. The spin polarization PM to AFM gives 
a considerable energy gain (~ 0.15 eV to ~ 0.3 eV per 
formula unit, depending on the magnetic phase), thus, 
polarization greatly lowers the energy regardless of the 
type of magnetic order. This is a sign that there is likely 
to be disordered local moment behavior above Tjv, al- 
though it is not so evident from the existing data.H Also, 
all AFM phases are more stable than the FM one. 

Now, we examine in more detail the competition 
among the different AFM phases. The AFM[ nl ] is worth- 
wile to be considered because it represents the usual kind 
of AFM order observed for NaCl-type structures. In- 
deed, a peculiarity of the AFM[ ln j phase (unless one 
considers noncollinear arrangements) is having all the 
six (fee) second neighbor moments antiparallel, a com- 
mon characteristic of most of NaCl .structures, as re- 
vealed by neutron diffraction studiesES. Nevertheless, 
we found the AFM[ ln ] clearly disfavoured with respect 
to the AFM^ 110 j and only slightly favoured on the FM 
phase. More interesting and somewhat surprising is the 
result for the AFM^ 110 j phase: we found it to be clearly 
the most stable of the cubic structures. It is character- 
ized by four spin-parallel and eight spin-antiparallel first 
neighbors, and all the six second neighbors have paral- 
lel spins (an arrangement very different from the other 
AFM phases here considered). But if the planar dis- 
tortion is added, the energetic order of AFM^ 110 j and 

AFM^ 110 j phases is reversed, and the AFM^ 110 j becomes 
the lowest, thus recovering the agreement with the exper- 
iment. Indeed, the distortion causes, for the AFM^ 110 j, 
the shortening of metal-metal distances between spin- 
antiparallel Cr atoms, that arguably causes an energy 
gain, and the stretching of spin-parallel Cr atoms in the 
(100) plane (producing an energy loss), resulting, as a 
matter of fact, in an appreciable net decrease of energy 
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(~ 0.07 eV/formula unit). The same cannot occur for the 
AFM^ 110 j , because no spin-antiparallel first neighbors are 
in the (100) plane, and thus the distortion does not allow 
any stabilization (actually we found an energy increase 
of ~ 0.01 eV/formula unit). 

Thus, the resulting physical picture is the following: 
if we consider the spin ordering in a constrained cu- 
bic symmetry, the AFM^q, phase, having the largest 
number of fee spin-antiparallel first-neighbors (8, against 
6 of both AFM| 110 j and AFM[m] phases) is favoured. 
But if we allow the symmetry to relax, the (100) pla- 
nar anisotropy of the AFM| 110 j phase causes some of the 
spin-antiparallel first neighbors to approach each other, 
producing a strong structural stabilization. We argue 
that in the cubic structure the bonds between spin- 
antiparallel metal atoms are under condition of large ten- 
sile stress (larger than that in the bonds between spin- 
parallel atoms, anyway) and the cubic to orthorombic 
transition relieves part of this stress. Then, the planar 
distortion plays a key role in the competition between 
AFM phases, and the stress relief is reasonably to be 
considered the ultimate driving force towards distortion. 

For magnetic moments (we show in the Table only val- 
ues referred to the experimental structures, being very 
close to that obtained for the theoretical ones) we ob- 
tain values ranging between m—1.9 \ib and 2.3 hb per 
Cr atom. In particular, for the AFM| 110 j phase our value 
rn=2.15 fiB is in fair agreement with the experimental 
m=2.36 fiB °f Ref. q, and the theoretical m=2.41 fis of 
Ref. ||. 

IV. ELECTRONIC PROPERTIES 

A. Effect of AFMf 110] Ordering 

The stabilization due to the magnetic ordering can be 
understood from the analysis of the density of states 
(DOS). Figure | shows the PM CrN DOS for the cu- 
bic structure in orthorhombic symmetry obtained from 
broadened eigenvalues from a regular mesh of k points in 
the irreducible zone, which is sufficient for our purposes. 
Unless otherwise indicated, in what follows all quantities 
refer to thfi^theoretical lattice constant. As pointed out 
previouslyEJ, among the series of 3d TMNs, CrN (as well 
as NaCl-type FeNE3) is distinguished for the high DOS 
at Ef. Our DOS agrees well with the full-potential lin- 
earized augmented plane wave calculation of Ref. [l2|. At 
least 95% of DOS at Ep comes from Cr d states, and the 
remaining contribution is due to the N p states. How- 
ever, the DOS peak at Ep is entirely due to Cr d states, 
indicating that the peak arises fcpm flat bands solely Cr 
in character. Singh and KleinEj have emphasized that 
this peak is Cr t<z g . As we show below, it is not possible 
to identify this flat band solely from symmetry line data. 

Figure [| shows total and partial DOSs of the cubic 
AFM| 110 j structure. The peak at Ep has disappeared: 



the AFM? 110 i ordering causes a large band splitting and 
the depletion (a pseudogap ~ 0.5 eV wide) in the DOS 
of a sharp region around Ep. 

Up- and down-polarized atoms are connected by trans- 
lation plus spin-flip symmetry operations, thus contribu- 
tions to DOS coming from the type of Cr and N atoms not 
shown in the Figure are equal to that shown but with up 
and down components interchanged. For N (middle panel 
of Figure |3|) magnetic effects are negligible, although not 
zero by symmetry: the magnetic moments are only ^0.01 
I^B- For Cr (bottom) the DOS around Ep clearly resem- 
bles the total one, hence the band splitting due to the 
polarization arises almost entirely from Cr d states. 



B. Fermi Surface Properties 

Notwithstanding the strong stabilization by moving 
states away from E^, AFM| 110 j ordering does not man- 
age to open a gap at the E^, and the system is still a 
(relatively low carrier density) metal, at variance with 
onea experimental characterization as a small gap semi- 
conductor but in agreement with the results of Browne 
et alS on very nearly stoichiometric samples. To provide 
a clearer picture of our prediction, in Figure ^ we show 
a band- resolv ed DOS including only the four bands (see 
Section £vd|) that cross or touch the Fermi level. These 
DOSs were calculated from a 528 k-points mesh, inter- 
polated using a Fourier spline techniqueca, and evaluated 
using the linear tetrahedron method. Full lines refer to 
the distorted (experimental) structure AFM| 110 j bands, 
while dashed lines refer to the cubic AFM^ 110 j phase. Ba- 
sically three bands contribute to the DOS at E^, and a 
fourth band has its minimum essentially at Ep. 

In Table |l| the band-by-band values of the band fill- 
ing, density of states at E^, and the Fermi velocities are 
given. From the totals we can evaluate the Drude plasma 
energies TiVL p : 



n 2 . 

p,j 



4TTe 2 N(E F )v 2 



(1) 



where j denotes any of the Cartesian axes. We obtain 
fi£l p ,x = 2.3 eV, Kl PtV — 2.0 eV, hil pz = 2.6 eV, so the 
non-cubic nature is evident but is only a ± 15% effect. 
This difference is due to the antiferromagnetism, not to 
the structural distortion whose effect we have found to be 
much smaller. These values of plasma energies are rather 
small for transition metal compounds, but are typical of 
metallic (not semiconducting or even semimetallic) be- 
havior. 



C. Analysis of Magnetic Structures 

We mention a technical detail about calculation of 
magnetic moments. In case of an AFM structure with 
aligned spins along high symmetry directions, it is very 
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suited for a plane-wave basis to extract magnetic mo- 
ments from the planar average of local magnetization, 
which is a more direct and simple procedure than the 
common way of projecting plane waves onto atomic-like 
orbitals. In the top panel of Figure ^| the planar aver- 
age of local magnetization along x for the FM (dashed 
line) and AFM| 110 j (full line) phase is shown. Peaks are 
in correspondence of atomic positions, and the integrals 
over interlayer distances give by construction the layer- 
by-layer magnetization, i.e. the magnetization per cou- 
ple chromium-nitrogen, or, in practice, the magnetization 
per Cr atom, being that of N atoms almost discardable. 
The idea to obtain contributions to some physical quan- 
tity coming from individual space regions, not accessible 
from its macroscopic (or 'integrated') value, using pla- 
nar and macroscopic averages has been sometime used. 
It leaded, for example, to introduce quantities such as 
energyEj and stressful density. In order to the integrals of 
the planar averages be well (uniquely) defined, a typical 
necessity is that the average is referred to an appropri- 
ate bulk symmetry plane (a critical case is, for instance, 
the average on (111) direction for a zincblebde structure: 
the atomic layers are not equally spaced and the volume 
attained to a single atom is undefined). We return at 
Figure || below. 

To complete the analysis of the DOS, in Figure || we 
present results for the FM phase of CrN in the exper- 
imental structure. As expected, the total DOS at the 
Ejr is dominated by the Cr d states. Ferromagnetism re- 
sults in the splitting of up and down d type bands, and a 
strong stabilization with respect to the PM phase occurs. 
The competition between antiferromagnetism and ferro- 
magnetism for metals is usually an intriguing question. 
In the present case, we can attribute the stronger AFM 
stability to a larger average spin splitting of Cr 3d states. 
It is clear from top panel of Figure || that for AFMj? 110 j 
there is larger magnetization on the Cr atoms. This trend 
should also be reflected in the band splitting. A measure 
of the average band splitting per atom is taken to be: 

(AV xc ) m = [ drAV xc (r) ■ m(r)/ f dr m(r) (2) 

i.e. the up and down exchange-correlation potential dif- 
ference AV XC , weighted by the local magnetization, inte- 
grated over the inter-planar volume Sj and normalized 
to the atomic moment. The integrand of Equation ^| is 
pictured in the bottom panel of Figure |^; the supralinear 
dependence of Kohn-Sham potential on magnetization 
leads to the difference between AFM^ 110 j (full line) and 
FM (dashed) AV xc (r) ■ m(r) in the regions around the 
atoms, resulting in an integrated value (over the inter- 
planar distance) of 2.58 eV for the distorted AFM^ 110 j, 

2.52 for the cubic AFM? 110] , and 2.02 eV for the distorted 
FM phase. The calculated splitting roughly reflects the 
magnitude visible in the Figures of DOS previously dis- 
cussed. In Figure U\ an expanded view of the partial DOS 



of Cr around E F is shown. Evidently, the AFM^ 110 j DOS 
has a minimun at Ep, whereas for the FM system E^ 
falls on the edge of the second peak of Cr d states that 
is unoccupied in the PM phase (see Figure ||). 

D. Band Shift from Magnetic Order and Lattice 
Distortion 

The band splitting mechanism is further clarified by 
analysis of band shifts. Figures || and ^|show the PM and 
the AFMj* 110 j band structures, respectively, calculated 
in the orthorhombic cell. The coordinates of the high- 
symmetry k-points are, in the orthorhombic reference 
system (see bottom panel of Figure ||), M=(7r/a,0,0), 
X=(0,n/b,0), L=(0,7r/b,7r/c), and W=hr/a,3-K/2b,0), 
where a, b, and c are given in Section O. In the cu- 
bic PM phase, at points M, X, W, and L all bands are 
doubly degenerate, due to the fee to orthorhombic band 
downfolding. Symmetry breaking due to the AFM^ 110 j 
ordering splits the degeneracy along the direction con- 
necting points M and W, whereas the lattice shear (a ^ 
90°) causes the band splitting at X and L. The former 
effect is the most effective in stabilizing the system. The 
splitting (~ 2 eV) strongly depletes the region around 
Ej?, leaving only few band crossings. It is particularly 
visible at M by comparing cubic PM and AFM? 110 , band 
structures. 

The effect of lattice shear produces a further stabi- 
lization, i.e. a further band opening, but is much less 
spectacular. It is emphasized in Figure ^[ band energies 
shown along M-T-X are those of the cubic PM phase, 
whereas along X-W-L the Figure reports values for the 
distorted structure, so that the splitting is visible at X. 
This effect is very small (its order of magnitude is only 
some tens of meV), and band structures of cubic and dis- 
torted AFM^ 110 j phase (Figures || and |To|) differ only in 
certain regions. 

An unusual feature that can be seen in Figure |t] is the 
coincidence (or nearly so) of E^ with four band crossings 
along the directions shown: two along the M-T line, one 
along T-X, and one at W. At general k-points away from 
these crossings, the bands will repel, hence they will not 
contribute to N(E^), but rather to a pseudogap at Ep- 

V. CONCLUSION 

Our first principles results for CrN confirm most of the 
available experimental data. The observed (110) double- 
sheet distorted structure is found to be the lowest in en- 
ergy, stable not only against paramagnetic and ferromag- 
netic phases, but also with respect to other antiferromag- 
netic arrangements. We found the system to be a weak 
metal, with large magnetic moments on Cr atoms. 

Although the fundamental features favouring a cer- 
tain spin arrangement with respect to another are still 
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to be understood, we are able to clarify the role played 
by the distortion by disentangling its effect and that of 
the magnetic ordering. Magnetic ordering is found to 
be dominant in determining both structural and mag- 
netic properties of the system, with respect to which the 
distortion can be seen as a small perturbation. Never- 
theless, distortion appears decisive to resolve the compe- 
tition among the different antiferromagnetic phases con- 
sidered here. Indeed, it causes a meaningful energy gain 
for the (110) double-sheet phase, possibly due to relieving 
tensile stress stored in the spin-antiparallel metal-metal 
bonds. 

Our calculations represent one of the first attempts 
to study magnetic systems within a framework of plane 
waves with ultrasoft pseudopotetials. All comparisons 
with previous available results confirm that the accu- 
racy of this methodology is comparable to that of the 
all-electron calculations. 
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TABLE I. Ef denotes energies (in eV/formula unit) for 
cubic structures calculated at the theoretical ao of the re- 
spective phases; Efs are relative to the experimental ao for 
all phases. Ef's and -E^'s are the same as before but now 
the structures are distorted by the observed change of angle 
a. All energies are referred to the that of the most stable 
structure, i.e. the AFM 2 110 ] phase with theoretical ao and 
distorted a. m's are magnetic moments (in //s/metal atom). 
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TABLE II. Band-by-band, and total values of the band 
filling (in electrons per Cr4N4 cell) DOS at the Fermi level 
(states/eV), and the rms values of the Cartesian components 
of the Fermi velocity (10 cm/s). The last line gives the total 
electrons and DOS, and the rms velocities at Ef for all four 
bands. 
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FIG. 1. Top panel: perspective view of the CrN or- 
thorhombic cell with antiferromagnetic spin arrangement con- 
sisting of alternating double ferromagnetic layers parallel to 
the [110] (following Corliss el al). Arrows denote the rela- 
tive spin directions on Cr atoms in the Neel state. Bottom 
panel: Top view of (001) plane. The dashed line indicates the 
orthorhombic cell with respect to the cubic structure. 
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FIG. 2. Density of states of cubic PM CrN. Top panel 
show the total DOS, the panels below show the major contri- 
butions resolved by angular momentum character. The sharp 
peak at E_f is entirely Cr d (note that the N p contribution is 
smooth through this region) . 
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FIG. 3. Density of states of cubic AFM CrN. Top panels 
show total and atom by atom DOS, middle and bottom panels 
some contributions resolved by angular-moment. 
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FIG. 4. Band-by-band contribution of the four bands 
around the Fermi level, for both cubic (dashed line) and dis- 
torted (full line) AFM^ 110 j structure. 
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FIG. 5. Top panel: the planar average of local magne- 
tization (in units /is/bohr) is pictured along the x = [110] 
direction of the orthorhombic cell for FM (dashed line) and 
AFM[ 110 ] (full line) phases. Bottom panel: planar averages 
(in Ryd'/is/bohr) of the corresponding exchange-correlation 
potential splittings (see text). 
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FIG. 6. Density of states of FM CrN calculated in the 
experimental structure. Top panel is the total DOS, lower 
panels show the individual contributions of each atom, as in- 
dicated in the Figure. 
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FIG. 7. Contributions of Cr to total DOS: comparison 
between AFM^ 110] and FM results. 
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FIG. 8. Top: band structure for cubic PM CrN in the or- 
thorhombic BZ. Energies between points M and X are for the 
cubic structure, where between X and L are for the distorted 
one, so that at X the minor band splitting due to the distor- 
tion is visible. Bottom: high symmetry points of the CrN BZ. 
Full line denotes the orthorhombic BZ, dashed line the cubic 
BZ and long-dashed line the fee truncated octhaedron. 
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FIG. 9. Band structure for the cubic AFMf 110] CrN. At 
variance with the PM structure, the double degeneracy at 
points M and W due to the cubic-to-orthorhombic downfold- 
ing is broken by the AFM^ 110 j symmetry, thus allowing the 
(mostly Cr d-like) bands to be shifted out from the Ef- 
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FIG. 10. Band structure for the distorted AFM? 110] CrN. 
Lifting of degeneracies can be observed, but there are no dra- 
matic changes with respect to the cubic AFM| 110 ] structure. 



